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Abstract 


For an understanding of a heat engine working in the microscopic scale, it is often necessary 
to estimate the amount of reversible work extracted by isothermal expansion of the quantum gas 
used as its working substance. We consider an engine with a movable wall, modeled as an infinite 
square well with a delta peak inside. By solving the resulting one-dimensional Schrodinger equation, 
we obtain the energy levels and the thermodynamic potentials. Our result shows how quantum 
tunneling degrades the engine by decreasing the amount of reversible work during the isothermal 
expansion. 
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I. INTRODUCTION 


A fundamental idea behind classical statistical mechanics is that we are living in a macro¬ 
scopic world, which means that we are unable to deal with microstates on a molecular level. 
This forces us to distinguish heat from work and to accept the increase of entropy on prob¬ 
abilistic ground, when we formulate the laws of thermodynamics. However, the border 
between our macroscopic world and the microscopic world of molecules is becoming vague 
due to technological developments over the last century. This has led to a natural question; 
i.e., what would it mean to thermodynamics if it became possible to access and manipulate 
microstates just as we do macrostates? The founders of statistical mechanics were already 
aware of this problem: For example. Maxwell imagined an intelligent being’s intervention 


Q. 


on a molecular level in his famous thought experiment, now known as Maxwell’s demon 
and expressed deep concern about the foundation of the second law of thermodynamics. 
The Szilard engine has been regarded as the simplest implementation of Maxwell’s de¬ 


mon with a single-particle gas j^. Its cycle consists of four processes; inserting a wall into 
the center of the box, measuring the position of the particle, expanding the quantum gas 
isothermally, and hnally removing the wall. Considering its microscopic nature, it is tempt- 


4 for other engines. In fact, there are more reasons than that: Zurek [^, when he discusses 
Jauch and Baron’s objection that the gas is compressed without the expenditure of energy 


ing to translate the cycle into quantum-mechanical language as has been done in Refs. and 


upon inserting the wall, argues that one has 


the Szilard engine jo]. According to Ref, 


:g 


:o resort to quantum mechanics to understand 
this ‘apparent inconsistency’ is removed by 


quantum-mechanical considerations. Therefore, in a sense, it is a matter of theoretical con¬ 
sistency. Kim and coworkers have presented a detailed account for its quantum-mechanical 
cycle in Refs. 7-^, with emphasis on the isothermal process. There seem to remain some 
subtle issues to settle, however, as shown in the debate on how to deal with the quantum 
tunneling effect through the wall in a three-boson case 10|. 

In this work, we describe the isothermal expansion of the quantum Szilard engine by 
considering a quantum gas conhned in a one-dimensional cylinder [ll, 12|. Although some 


previous studies, such as Refs. 


13-15 


, assume that the wall is impenetrable, we will relax 


that assumption because it is, strictly speaking, experimentally infeasible. Starting from the 
Schrodinger equation, we calculate the energy levels and thereby obtain how much work can 
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be extracted by an isothermal reversible expansion. We find the following: T 


reversible work is significantly rednced compared to the claims in Refs. 0 and [7 
full thermal eqnilibrinm inside the box, as is usnal in describing isothermal expansion. This 
conclusion should not be affected by either the wall’s insertion or the measurement before the 
isothermal expansion because thermal equilibrium does not remember its history. Moreover, 
if we put two or three fermions in the box, the free-energy landscape is not monotonic with 
respect to the wall’s position, so that we should sometimes perform work to expand the gas. 
These effects have not been explicitly discussed in previous studies. 

This work is organized as follows: An explanation of our basic setting in terms of the 
Schrodinger equation is given in Section [ITl A numeric calculation of the energy levels and 
the free energy and a comparison of a single-particle case to the two- and three-particle cases 
is presented in in Section IIIII This is followed by discussion of results and conclusions. 


le amount of 
if we assume 


II. PARTICLE IN A BOX 


Consider the following one-dimensional potential landscape with size L: 


V{x) 


a6{x — pL) a 0 < X < L 
oo elsewhere. 


( 1 ) 


where a is the strength of the delta potential and p G (0,1) specifies its location 
expect the Schrodinger equation 



. We 


2m dx'^ 


V{x) 


'0(x) = Efj^x) 


( 2 ) 


to have a solution with a sinusoidal form. The boundary conditions at x = 0 and x = L are 
satisfied if we set 


fj^x) 


A sin kx if 0 < x < pL 

B sin k{L — x) if pL < x < L, 


( 3 ) 


where the coefficients A and B, as well as the wavenumber fc, are assumed to be nonzero. 
The eigenfunction i/^(x) has the following properties: First, it is continuous over [0,L]; i.e.. 


A sin pkL = B sin qkL 


( 4 ) 
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with q = 1—p. Second, the change in the derivative of with respect to x around x = pL 
is obtained by using 


ppL+e 

I pL—e 

from which it follows that 


lim 

e^O 




, 2m 
dx + — 


ppL+e 


E%lj{x)dx 


ppL+e 


a'ip{x)S{x — pL)dx 


_J pL—e 


/ pL—e 


= 0. (8) 



pL+e 


difj 

dx 


pL—e. 


2ma 




( 6 ) 


here and multiplying both sides by sin qkL to use Eq. (jl]), we obtain 


—k sin kL = sin pkL sin qkL. 


(7) 


This formula has also been derived in Ref. Il2| in a different way. Equation o can be 
numerically solved by using the Newton-Raphson method to yield the allowed values of k. 
The method works as follows: First, let us dehne f{k) = k sin kL + sinpkL sin qkL, 

and look for its zeros to solve Eq. ([7]). Starting from fc = 0, we check whether the sign of 
f{k) changes with increasing /c by a sufficiently small amount, say, dk <C 1. Every time the 
sign changes between k = k* > 0 and k* + dk, we run the following iteration starting from 

fcW = k*: 


kiJ+L 


kU) _ 


/[fcW] 

f'[k(Ly 


( 8 ) 


where j = 0,1, 2,.... When k^^^ has converged to a stationary value k^°°\ i.e., if k^^'^ = 
within numerical accuracy, we take it as a solution of Eq. ([7]) and repeat the procedure 
by checking the sign of f{k) from k = k^°°'^ + dk. Each of these solutions is identihed 
with the wavenumber kn of the nth eigenmode in an ascending order. For the application 
of the Newton-Raphson method numerically, it is convenient to express the quantities in 
dimensionless units, so we divide Eq. (|2]) by Eg^'^ = h^7r^/(2mL^), the ground state energy 
for a = 0, to obtain 


A 

dp'^ 


^{p) + ^{p) 


V’(h) = 


where p = nx/L. Then, the potential inside the box is expressed as 
^(p) = -pL) = a6{p-p7i). 


(9) 


( 10 ) 
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FIG. 1: (Color online) (a) Lowest wavenumbers in units of vr/L as a function of p, and (b) the 
corresponding energy levels in units of /(2mL‘^). The other two panels show the 

probability densities for the three lowest eigenmodes (c) when the wall is in the middle, p = 0.5, 
and (d) when it has moved to p = 0.4. The strength of the potential is set as a = 10. 

with a = 2amL/{fi?n) = a{L/n)/Ef\ We can also define a dimensionless wavenumber 
K = Lk/ir to rewrite Eqs. (I7j) as Ksin ktt + asinp/tTr sin = 0. A little algebra shows that 
the energy level is given by 

( 11 ) 

where the second term represents the potential energy coming from the overlap between the 
wavefunction and the delta peak potential. The numerical solutions for these k and e are 
depicted in Fig. [H together with the probability density plots. 


, 2a 
e = -\ - 

71 


CSC pKTT p 


2k7T 


sin 2pK7T + CSC qK7T q 


2K7r 


sin 2qK,7i 


III. REVERSIBLE WORK 

A. Single Particle 

We are interested in how much work can be extracted by performing a reversible isother¬ 
mal process with this system. Let e„ denote the nth energy level obtained by using the 
Newton-Raphson method with n = 1,..., M so that the ground-state energy is denoted by 
Ci. For a given value of temperature T = where ks is the Boltzmann constant, we 
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FIG. 2: (Color online) Change in the free energy F in units of = h^7r^/(2mL^) when the wall, 
starting from the center of the box, moves to a new position p. The height of the wall is fixed at 
a = 10, and each curve represents a different T in units of Eg^'^/ks- The number of energy levels 
considered is kept as M = 10^ throughout the free-energy calculations. 


run the calculation to get /3eM S> 1. The partition function for a single quantum particle 
can then be obtained as 

M 

Z= lim (12) 

M^oo ^ 
n=\ 

The free energy is F = —ksTlnZ, accompanied by a differential form dF = —SdT — 
PdV, where S', P, and V are the entropy, pressure, and volume, respectively. The classical 
concepts, such as pressure and work, need careful consideration in quantum mechanics [l^ . 
To calculate work, we begin with the occupation probability of the nth level given by 


Pn 


V e~P‘ 

Z-^n 


-ksTidZ/den) 

z 


A 

den 


{-keT In Z) 


dF 

den 


(13) 


The change in the internal energy U = is then expressed as dU = endpn + 

Pnden = dQ + dlT, where dQ and dW are differentials of heat and work, respectively. The 
quantum thermodynamic work is identihed with dW = because heat is usually 

involved with changes in occupation probabilities for given energy levels whereas work is 
performed when the energy levels themselves change [71] . Consequently, the amount of 
reversible work during this process is W = J ~ I ~ which is consistent 


with the classical case. 

Figure [2] shows AP resulting from moving the wall, which was initially located at the 
center, p = 0.5. Suppose that T is so low that T < /ks- When the wall is at p = 0, 
only the ground state contributes to the summation in Eq. flT^ . which means that Z{p = 
0) \ When the wall is at p = 0.5, the ground-state energy equals 4Pg°\ and the 
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first excited state also contributes to the partition sum in Eq. flT^ because it lies close 
to the ground state [see Fig. [I][b)]. The partition function is, therefore, approximated as 
Z{p = 0.5) ~ \ The free-energy difference in this low-T region is, thus, given 

as T(0) — F(0.5) ~ —3E® + ksT\n2. When T = /ks, for example, the free-energy 
difference roughly amounts to —2.3Eg^\ which explains the size of the free-energy drop in 

Fig. H 

As T increases, a plateau develops near the center of the box because E does not respond 
much to the wall’s position. However, the drop in E at the end of the process increases with 
increasing T if T ^ 0(1) in units of E^'^/ks. Our question is how it grows with T. Let us 
take a —)■ oo and 7 = /3h^7r‘^/{2mL‘^) <C 1. Then, it is the partition function for the initial 
and hnal wall positions can be approximated. The latter case of the hnal wall position is 
estimated as 


Z(p = 0) = Y.' 


,- 7 ^ 


n=l 

poo 


'0 


' I 

exp(- 7 n^)dn- 


g-7r),2 _ ^-7(n+l)2 


n=0 


TT 1 

47 ~ 2’ 


( 14 ) 

( 16 ) 

( 16 ) 


where the summation in Eq. (1T5|1 is a hrst-order correction to the integral approximation. 
Likewise, the former case of the initial wall position gives 


Z{p = 0.5) = 

n=l 


2 



exp (— 4771 ^) 0 ? 




(17) 


where the factor of 2 in front of the summation is due to the fact that every pair of adjacent 
energy levels becomes degenerate in the limit of a —>■ cx). This effect cancels the factor 
of 4 inside the exponential arising from the reduced volume L/2. Although AE does not 
completely vanish due to the correction, the important point is that Z(0.5)/Z(0)_^for 7 <C 1 
is far less than 2, in contrast with the claim that AE = —kBT\yi2 in Refs. l5| andlT]. Because 
a must be hnite in any experimental situation, the particle can, in full thermal equilibrium. 


extracted work. On the 
5|, I?! have assumed that 


be observed on either side of the box, which reduces the amount o: 
other hand, the previous studies on the quantum Szilard engine 
the process is performed within a shorter time scale than required for tunneling through the 
wall and, therefore, concluded that Z{p = 0.5) ~ f exp{—4''yn‘^)dn without the factor of 2 in 


6 














p 


p 



p 


LO 

0 


0 

(d/y 





T=1 - 

1 

-5 


3 

5 


(i^^'Three 

10 . 

Fermion^O . 

Uj 


0 0.1 

0.2 0.3 0.4 0 


p 


FIG. 3: (Color online) Free-energy changes in the presence of multiple non-interacting quantum 
particles: (a) bosonic and (b) fermionic cases with two particles. With three particles, (c) the 
bosonic case remains qualitatively similar, but (d) the fermions behave differently due to the Pauli 
exclusion principle. 

front. Strictly speaking, their engine nndergoes an isothermal process only when partially 
eqnilibrated to maintain the conhnement and extract W = ksT In 2. 


B. Two- and Three-particle Cases 


If two qnantnm particles are in the box, we should consider their symmetry, i.e., whether 
they are fermions or bosons. Without any consideration of the spin degree of freedom, the 
elements of the density matrix for the two particles are expressed in the bracket notation as 


follows: (e'^e 2 


,-/3H 


6162 ) = e ^hi+^ 2 ) ± where 6 is the Kronecker delta 


symbol, and we have a plus sign for bosons and a minus sign for fermions. The two-particle 
partition function is then obtained by taking the trace operation: 

1 


Z = Ti' 1 e-'’") = ^ (1 ± 4.„). 


(18) 


«=:' P 


,-/SH 


The same applies to the three-particle case. The density-matrix elements are written as 

), and the partition function reads 

(1 i *^eie 2 ^ ^£2^3 ^ *^£ 3^1 T 2*^eie2*^e2e3) ' 


616263) = e / 3 (^i+" 2 +" 3 ) ± 


Z = Tr = g-/3hl+e2-l-e3) 

£ 1 ,^ 2 ,es 
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Figure [3] shows the free-energy differences for the two- and the three-particle cases. The 
bosonic cases are qualitatively similar to the single-particle case. However, we see very 
different behavior for the fermionic cases where the Pauli exclusion principle is in action: 
At low T, we need to perform positive work to move the wall from the center to a certain 
position p [Figs. El^b) and (d)]. Obviously, the reason is that it costs free energy to place 
two fermions close to each other. 

IV. SUMMARY AND OUTLOOK 

In summary, we have considered an isothermal expansion process of a quantum gas, 
taking the tunneling effect into consideration. We have found that the amount of reversible 
work is smaller than In 2 in the high-T region if the system is fully equilibrated at every 

moment of the isothermal expansion. The difference from Refs. 0 and Q arises because they 
have separated the time scale of tunneling from that of the partial equilibration that occurs 
on only one side of the wall. Although it was not explicitly stated in previous studies, 
this separation may be a plausible assumption for the following reason: As we increase 
the potential height a, we may well expect the time scale for tunneling to grow whereas 
the partial equilibration before tunneling is achieved within a finite amount of time. The 
quantum Szilard engine will show its expected performance only between these two time 
scales. The question is, then, how large a value of a one should have to ensure the separation 
of the time scales, which will be pursued in our future studies. 
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